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Abstract. We construct a nonequilibrium theory for the charge transfer through 
a diffusive array of alternating normal (N) and superconducting (S) islands com- 
prising an SNSNS junction, with the size of the central S-island being smaller than 
the energy relaxation length. We demonstrate that in the nonequilibrium regime 
the central island acts as Andreev retransmitter with the Andreev conversions at 
both NS interfaces of the central island correlated via over-the-gap transmission 
and Andreev reflection. This results in a synchronized transmission at certain reso- 
nant voltages which can be experimentally observed as a sequence of spikes in the 
differential conductivity. 



1.1 Introduction 

An array of alternating superconductor (S) - normal metal (N) islands is a fun- 
damental laboratory representing a wealth of physical systems ranging from 
Josephson junction networks and layered high temperature superconductors 
to disordered superconducting films in the vicinity of the superconductor- 
insulator transition. Electronic transport in these systems is mediated by An- 
dreev conversion of a supercurrent into a current of quasiparticles and vice 
versa at interfaces between the superconducting and normal regions §. A fas- 
cinating phenomenon benchmarking this mechanism is the enhancement of the 
conductivity observed in a single SNS junction at matching voltages constitut- 
ing an integer (m) fraction of the superconducting gap, V — lAjiem) P-pd| 
due to the effect of multiple Andreev reflection (MAR) [p^[-p| . The current- 
voltage characteristics of diffusive SNS junctions were discussed in great detail 



2 N. M. Chtchelkatchev^'^'^ T. I. Baturina^-'*, A. Glatz^ and V. M. Vinokur^* 

M = 1 : cVm = 2A 




n = 2: ey,o, = A 




Fig. 1.1. Diagrams of tlie SAT processes for tlie first, n — 1, (a)-(b), and second, 
n = 2 (c)-(f), subliarmonics of tlie resonant singularities in dl/dV described by 
Eq. ( [i.l| ) for tlie SNSNS junction with the normal resistances ratio R1/R2 — 3/4 
(depicted through 3/4 ratio of the respective lengths of the normal regions). The 
thick solid lines represent the quasiparticle paths starting and/or ending at the 
points of singularity in the density of states at energies e — itA at the electrodes 
Sl and Sr. The dashed solid lines show paths starting and/or ending at the edges 
of the gap of the central island Sc. The circle denotes the over-the-gap Andreev 
reflections at the electrodes. The paths ABCD [panel (a)] and D'C'B'A' [panel (b)] 
correspond to the electron- and hole trajectories, respectively. Synchronization of 
the energies of the incident and emitted quasiparticles at points B and C (B' and 
C) is shown by arrows. SAT is realized by trajectories passing through the singular 
points e = ±A of the central island Sc and including over-the-gap transmissions 
and Andreev reflections. Trajectories synchronizing other transmissions across Sc 
and those of higher orders are not shown. Note, that voltage drops eVi = 6A/7 and 
eV2 = 8A/7 [panels (a)-(b)] {eVi = 3A/7 and eV2 = 4zi/7 [panels (c)-(f)]) are not 
MAR matching voltages of individual SlNiSc or ScN2Sr parts. 
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in Rcfs. p5|,pi|. Further developments were obliged to studies of large arrays 
comprised of many SNS junctions p7|-pl[|. Experimental results, especially 
those obtained on the multiconnected arrays ||l^,[8[|2^, indicated clearly that 
singularities in transport characteristics cannot be explained by MARs at in- 
dividual SNS junctions only and that there is evidently a certain coherence 
of the Andreev processes that occur at different NS interfaces. These findings 
call for a comprehensive theory of transport in large SNS arrays. 

In this article we develop a nonequilibrium theory of electronic transport in 
a series of two diffusive SNS junctions, i.e. an SNSNS junction and derive the 
corresponding current-voltage characteristics. We demonstrate that splitting 
the normal part of the SNS junction into two normal islands that have, in gen- 
eral, different resistances and are coupled via a small superconducting granule, 
Sc leads to the nontrivial physics and emergence of a new distinct resonant 
mechanism for the current transfer: the Synchronized Andreev Transmission 
(SAT). The main component of our consideration is a nonequilibrium circuit 
theory of the charge transfer across Sc. [The symmetric case with the equal 
resistances of the normal parts was discussed in detail in Q . Unfortunately 
the technique developed there does not allow straightforward generalization 
onto a nonsymmetric case.] 

In the SAT regime the processes of Andreev conversion at the boundaries 
of the central superconducting island are correlated: as a quasiparticle with 
the energy e hits one NSc interface, a quasiparticle with the same energy 
emerges from the other SqN interface and enters the bulk of the normal is- 
land (and vice versa, see Fig. 1). This energy synchronization is achieved via 
over-the-gap Andreev processes , which couple MARs occurring within the 
each of the normal islands and make the quasiparticle distribution at the cen- 
tral island essentially nonequilibrium. Effectiveness of the synchronization is 
controlled by the value of the energy relaxation lengths of both, the quasipar- 
ticles crossing Sc with energies above A, and the quasiparticles experiencing 
MAR in the normal parts. The SAT processes result in spikes in the differ- 
ential conductivity of the SNSNS circuit, which appear at resonant values of 
the total applied voltage Vtot defined by the condition 

2A 

14ot = — (1.1) 

en 

with integer n, irrespectively of the details of the distribution of the partial 
voltages at the two normal islands. 

The article is organized as follows. In the Section 1.2 we define the system, 
a diffusive SNSNS junction which will be a subject of our study. Section 1.3 is 
devoted to introduction and description of the employed theoretical tools: the 
electronic transport of the system in the resistive state is given by the Larkin- 
Ovchinnikov equation in a form of matrix equations for the Green's functions 
taken in Keldysh representation. In Sections 1.4-1.8 we construct an equivalent 
circuit theory for an SNSNS junction resulting in the recurrent relations for 
the spectral current flow in the energy space. In Section 1.9 we present the 
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original numerical method enabling us to solve the recurrent relations for the 
spectral current and obtain the I-V characteristics for the SNSNS junction. 
The obtained results are discussed in Section 1.10, where we demonstrate, in 
particular, that the SAT-induced features become dominant in large arrays 
consisting of many SNS junctions. 



1.2 The system 

We consider charge transfer across an SlNiScN2Sr junction, where Sl, Sc, 
and Sr are mesoscopic superconductors with the identical gap A; the 'edge' 
superconducting granules Sl and Sr play the role of electrodes, and Sc is the 
central island separating the two normal parts with, in general, different nor- 
mal resistances. We discuss the common experimental situation of a diffusive 
regime where the most of the energy scales are smaller than h/r, where r is 
the impurity scattering time. We assume the size, Lc, of the central island 
to be much larger than the superconducting coherence length ^, hence pro- 
cesses of subgap elastic cotunneling and/or direct Andreev tunneling [ p3| do 
not contribute much to the charge transfer. In general, this condition ensures 
that Lc is large enough so that charges do not accumulate in the central 
island and Coulomb blockade effects are irrelevant for the quasiparticle trans- 
port. At the same time Lc is assumed to be less than the charge imbalance 
length, such that we can neglect the coordinate dependence of the quasipar- 
ticle distribution functions across the island Sc. Additionally, the condition 
ig ^ Lc, where is the energy relaxation length, implies that quasiparticles 
with energies e > A traverse the central superconducting island Sc without 
any noticeable loss of energy. The normal parts Ni and N2 are the diffusive 
normal metals of length Li 2 > and L12 > L^, L^ — ^ KD^/e, where 
is the diffusion coefficient in the normal metal. We assume the Thouless 
energy, i?Th = ^^n/-^i,2j to be small, i^xh ^ A, and not to exceed the char- 
acteristic voltage drops, iSxh < eVi^2- If these conditions that define the so 
called incoherent regime ||l^ are met, the Josephson coupling between the su- 
perconducting islands is suppressed. And, finally, we let the energy relaxation 
length in the normal parts Ni and N2 be much larger than their sizes, so that 
quasiparticles may experience many incoherent Andreev refiections inside the 
normal regions. 



1.3 Theoretical formalism 



The current transfer across the SNSNS junction is described by quasiclassical 
Larkin-Ovchinnikov (LO) equations for the dirty limit |^ , |25| : 

-j[ifeffO,G] = VJ, J.n=-i-[G3,GJ, (1.2) 

za^H 
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where HcS = ^ii^zdt — f<7o + J — DG o VG is the matrix current, the 
subscripts "S" and "N" denote the superconducting and normal materials, 
respectively, "o" stands for the time-convolution, ai (i = x,y,z) are the Pauli 
matrices, operating in the Nambu space of 2 x 2 matrices denoted by 'hats', 
A = i(7x Im A + i&y Re A, and R is the resistance of an NS interface. The 
diffusion coefficient D assumes the value D,^ in the normal metal and the 
value Dg in the superconductor, and f is the electrical potential which we 
calculate self-consistently. The unit vector n is normal to the NS interface 
and is assumed to be directed from N to S. The momentum averaged Green's 
functions G{r, t,t') are 2x2 supermatrices in a Keldysh space. Each element 
of the Keldysh matrix, labelled with a hat sign, is, in its turn, a 2 x 2 matrix 
in the electron-hole space: 

r is the spatial position, t and t' are the two time arguments. The Keldysh 
component of the Green's function is parametrized as |Q : ~ G^ o / — / o 
G^, where / is the distribution function matrix, diagonal in Nambu space, 
/ = diag [1 — 2no, 1 — 2nh]j '^o(h) is the electron (hole) distribution function. 
In equilibrium ne(h) becomes the Fermi function. And, finally, the Green's 
function satisfies the normalization condition G^ = 1. 

The edge conditions closing Eqs. (^]^) are given by the expressions for the 
Green's functions in the bulk of the left (L) and right (R) superconducting 
leads: 

GL(R)(t,i') = e-^^L(R)tf3AG^(i _i')e'^L(R)t'f3//i^ 

the chemical potentials are — and fi^ = eV. Here, Go(<) is the equilib- 
rium bulk BCS Green's function. 

The current density is expressed through the Keldysh component of J as 

I{t,r) = !I^TraJK(i,t;r) = ^J de[Ue)+h{e)] , (1.4) 

where the spectral currents Ic and /h representing the electron and hole quasi- 
particle currents, respectively, are the time Wigner-transforms of top- and 
bottom diagonal elements of the matrix current j'^'. 

On the normal side of the superconductor-normal metal interface, the 



Keldysh component of Eqs. (1.2) yield the conservation conditions: 



V/e(h) = 0, (1.5) 

/o(£) = {Dp{e + u)Vne(£) - An(£ + u)Vnh(£ + 2u)} , (1.6) 
/h(£) - {Dpie - u) Vuhie) - A„(£ - u)Vn^{e - 2u)} , (1.7) 

where u is the electrical potential of the adjacent superconductor, -Dp(m) = 
iD^±D+)/2, 



6 



N.M. Chtchelkatchev^'2'^, T. I. Baturina^-"*, A. Glatz^, and V. M. Vinokur^* 




Fig. 1.2. Effective circuit for tlie boundary between the normal metal and the 
superconductor. Kirchhoff laws where the role of the potential in the nodes is taken 
by the electron- and hole distribution functions give the boundary conditions for 
the LO equations, (a) A general nonequilibrium case, (b) Equivalent circuit for an 
equilibrium case where quasiparticle distribution functions in the superconductor 
are the Fermi-functions, np, then ne{e) = nF(e -|- m) = nh_{s + 2u), used in The 
superconductor electrical potential at the boundary is equal to u. 

4 (1.8) 
D_{e) = -Tr[i-cr"G='(e)a" ^^(e)] , 

and the trace is taken over components in the Nambu-space. In the bulk of 
a normal metal, G'"''^'(e) — ?► ±&z and Z)+ w £>_ w 1, so Ic = CTNVric and 



1.4 Circuit representation of the boundary conditions 

We start the construction of the circuit theory with the corresponding for- 
mulation of the boundary conditions for the distribution functions at the in- 
terface between the normal parts and the central superconducting island. We 
consider a stationary situation where the applied voltage does not depend on 
time. Then the Green's functions can be parameterized near an NS interface 
as follows: 

[G^]j(e,e') =az^e-e' cosh9j{e) + 

&+ Se-£'+2u sinh6'j(e) - &- Se-£'-2u sinh6'j(e) , (1.9) 
G* = -a,(G-)ta,, 

where a± = CTx ± icTy, j = S, N, and u is the electrochemical potential. 
The effective diffusion coefficients are correspondingly £)+ = cos^ Im 9 and 



£>_ — cosh^Re6'. When deriving Eq.(L9), we have used the condition that 



the Josephson coupling between the superconducting islands in the junction 



is suppressed. The proximity effect results in an additional term in Eq.(1.9) 
proportional to d{e — e' — 2{u' — u)), where u' is the potential of the adjacent 
superconductor involved. 
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Taking the Keldysh component of the boundary term in Eq. ( p_.2| ) we derive 
the boundary conditions for the currents /o(h) at the NS interface, which 
assume the form of Kirchhoff's laws for the circuit shown in Fig. ^]^(a). The 
electron and hole distribution functions take the role of voltages at the nodes. 
The equation for an electronic spectral current flowing into the lower left 
corner node: 



p(lC) 



(s + u) 



iif'(£ + 2u)-<^'(£) 
[-i?<,^°'(e + u)] 



<"(e + 2u)-n(i'(£) 



(1.10) 



The equation for the hole current going into the to p lef t node of the circuit 
of the Fig. ^]^(a) is easily obtained analogously to ( I.IC ) with the aid of the 
additional transformation e — >■ e — 2u i.e. by shifting all the energies over —2u: 



r <"'(£)-<'(£) ^ 2^) -<(£) , 



n(iC) 



n<^'(e-2u)-ni\''(e) 



(1.11) 



In an equilibrium the quasiparticles in the superconductor follow the Fermi 
distribution, then n<,'^'(e) = n^is + u) = n{f' (e + 2u) and Eqs.( 1.10 )-( 1.11 ) 
reduce to 



np (e - u) - n<i) (e) n'^^ [e + 2u) - n<^> (e) 



(e + u) 



(1.12) 



where the interjacent resistances are defined as R^l^){e) = {RZ^{e) ± 



R-\e)}/2; here l/i?±(£) = [N2 Ni T M^]/ R, Nj{e) = Recosh6'j, 
— sinOj and j — 1,2 labels the different sides of the in- 
terface. 



The circuit representation of Eq.( 1.12 ) is shown in Fig. L2(b). This is the 
boundary conditions and the corresponding circuit used in Ref. |]l5| . 



1.5 Conductance renormalization procedure 



We consider the normal metal between the left superconducting lead Sl 
and the superconducting island, Sc, see Fig. 1.3. The boundary conditions, 
Eqs.(l.lC)-(l.ll), relate electron and hole distribution functions at the right 
NS and left NS interfaces. Below we relate electron and hole distribution func- 
tions at a; = and x — d building the effective circuit, where d is the length 
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Fig. 1.3. Illustration of the spectral currents flow in the normal metal [white area] 
surrounded by the superconductors [grey area]. The black boxes at the interfaces 
encode the boundary conditions picture like it is shown in Fig.[L.2|. 



of the normal layer. At the first step we neglect the proximity effect change 
of the junction resistance and take D± — 1 everywhere in the normal layer. 
Then, — Vn(,(h), An^jh) — and therefore 

4(h) (e) = [nc{h){d,e) - ne(h)(a; = 0,e)]/i?i , 



(1.13) 



where i?i = ai/dis the normal resistance of the Ni-layer. Eq.( 1.13 ) resembles 
the Ohm law for the resistor i?i , but where the role of the v oltages play the 
distribution functions at the ends of the resistor. Eq.(1.13) is approximate 
because it neglects the proximity renormalization of the normal layer conduc- 
tivity It was shown in Ref. [|l^ for a SNS junction that the replacement of 
ne(h)(a^ = {0, d}) by the properly chosen proximity renormalized distribution 
functions makes Eq.(1.13) accurate. We show below that this idea is appli- 
cable when electron and hole distribution functions in the superconductors 
essentially deviate from the Fermi functions and when the electron and hole 
currents can not be in general related by a shift of the energy like in SNS 
junction. 

At the left NS-interface the spectral currents, /c(£) and /h(e + 2u) are 



related by the Andreev process, see Fig. 1.2a. It follows from Eqs.(1.5)-(1.7) 



that the combination of the quasiparticle currents, /j!°'(e) = /o(e) ± /h(£ -t- 
2u) = CTi-D^'^' (e + u)Vn^'^'(e -I- u), conserve in the normal metal: V/j!°' = 0. 
Integrating the last equation over x we get, 

dx 



7-(lC) 



= CTi [n±{d) - n±{x)], 



(1.14) 



where ni^*^'(e) = nc{e) ± nh(e + 2u). Eq.( |l.l4| ) can be equivalently rewritten: 



J-(IC) 



(d-x) 
\d)=n[ 



ox \n''^\d) 



,(1C) 



{x)l 



\d) 



mi^°'/r'(e + u), 



(1.15) 
(1.16) 



where 
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Here the variable x occupies the domain ^ 2; ^ ^ where the Cooper 
pair wave functions from the left and right superconductors, see Fig. 1.3, do 
not overlap. At these values of x, the angle 9{x) — >■ 0, D± {x) — » 1 and we can 
therefore substitute a: by in the integral written in Eq.(1.14). 

Taking into account that I±'^\s + u) = Ic{s) ± /h(e + 2m) we finally get 
the following important result: 



{d~x)I,{s) = <J,[n^^-\d)-n^^-\x)], 



(1.18) 



where 



(d) = n<^^) (d) - /e(e) m<^^' (e + u) - his + 2u) m\^^> (e + u) . (1.19) 



Here m^'^^'^ = [n 

Applying the procedure, Eqs.(1.14)-(1.19), to the left NS interface in 
Fig. 1.3, we find for ^ ^ — ^n- 



a:/e(£) -ai[<^'(:r)-n<^^'(a: = 0)] 



(1.20) 



where 



ni'^\x = 0) = n<^^'(.T = 0) + /e(£) m<^^'(£ + u') + 

h{e + 2u')m\^'->{e + u'). (1.21) 



c(h) 



1 dx. 



(1.22) 



Eqs.(|1.18|)-( pr20|) sh ow how nP depends on X in the central part of the 
normal layer in Fig 1.3. Eq.(1.18) must be consistent with Eci.(1.20). The only 
way to satisfy this condition is the following one: 

n^^'^>{d)-n'^''^Hx = 0) 



Ri 



(1.23) 



where i?i = d/ai is the normal resistance of the N-layer and we used that 
n^^^^x) = n'g'^'>{x). The condition Eq.( 1.23| ) resembles the Ohm law. It allows 
to relate the distribution functions at x = and x — d. 
Similar condition holds for h: 



Ri 



(1.24) 
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where 



<°'(d) - m<^^'(e - u) h{e - 2u) - m<^°'(e - u) h{e) , (1.25) 



^<^^^' = n'^^' + mj^^^'(e - u') I,{e ~ 2u') + m<^^'(£ - u') h{e)] 



(1.26) 



The last step is the formulation of the boundary conditions at the NS inter- 
faces in terms of the distribution functions with bars. Using the I± notations 
we can rewrite the boundary conditions, Eqs.(Ll^-(l.ll), in the compact 
form 



(e + u) 



(1.27) 



Then it follows from Eg. (1.15) that we can write: 



(C) -(IC) 

n\_' — n_L. 



RPis)=m'l^\e)+R'P{e). 



(1.28) 
(1.29) 



The same form has the boundary condition at x = 0: 



R'P{e)^mP{e) + R'-l^\e). 



(1.30) 
(1.31) 



It follows that the physical meaning of m± terms is the proximity effect con- 
tribution to the NS interface resistance, see p^ . 

It is more convenient to work with the boundary conditions for /o(h) rather 
then with those for I±. Then one can use Eqs.( 1.10 ),( 1.11) but with n'-^^^ 



(IL) 

e(h) 

R^] 



"e(h) and i?Q(p) 



^Q(p)' ■^hfire, for example, -Rq(p) 



(IL) 

e(h) 

2^(1C)^(1C)/J^(1C)^ 



1.6 Retarded and advanced Greens functions evolution 
in normal metals and superconductors 

Having formulated the boundary conditions for the distribution functions we 
turn now to advanced and retarded Greens functions behaviors. 

Normal layers in experimental SNS junctions and SNS a rray s, see Ref. 
T7[|l9|, connect with superconductors like it is shown in Figs. |l.4| . The junc- 



tions of this type are usually referred to as "weak-links". [E7, 28] Boundary 
conditions for retarded and advanced Greens functions, Eq.(|l.2), can be sim- 
plified in this case: retarded and advanced Greens functions at superconduct- 
ing sides of NS boundaries can be substituted by Greens function from the 
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bulk of the superconductors. These "rigid" boundary conditions approxima- 
tion is reasonable because the magnitude of the current is much smaller than 
the critical current of the superconductor [this is assumed] and the current 
entering the superconductor from narrow normal metal wire with the width 
comparable with the Cooper pair size. There are also other cases when the 
rigid boundary conditions are correct, for example, if the NS boundary has 
the small transparency due to, e.g., an insulator layer at the NS interface. 

The recipe telling how one should evaluate d{x) in the normal metal near 
the NS boundary, where the rigid boundary conditions hold, can be taken 
from e.g. Ref. |1^, and we reproduce briefly their result for the c omp leteness. 
We will write down 0{x) near the right NS boundary (see Fig.^^s] ) taking 
u' = 0. In the superconductor, 9s = atanh(zi/£), where A is the gap. The 
value of 6'n = 9{x = 0) in the normal metal side should be found from the 
equation: 



wJ sinh(0^-gs) + 2sinh% ^0, (1.32) 

y e + i/2Tcr 2 



where W — Ra/Rns, Ra = Ca/^n, (Ca — \J D j A) is the resistance of the 
normal metal layer with the width ^a. Here is the pair breaking rate [ p9[ 
[e.g., induced by electron-phonon or electron-electron interactions] and i?Ns is 
the normal resistance of the interface. Then the solution for 9{x > 0) is the 
following: 

tanh — = exp I p|tanh— , (1.33) 

4 V ^.ViJ 4 



where £,£ — y/D/[2{e + i/2Ta)]. Th e effe ct ive c onductances g± should be ex- 
pressed through 9 found from Eqs.(1.32)-([l.33|). 




Fig. 1.4. Typical experimental array of SNS junctions [p^p^. This type of the 
link enables us to use the rigid boundary conditions for the retarded and advanced 
Greens functions. 
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l„{e+2u) 



(a) ni'\8+2u) R,$'^'(6+u) 


nj'^'(£+2u) 


<"{e+u) ni"(£+2u) 


Ih(^+2u) 

+ A 

G 




Sc 




Ih(e+2u) 

o 


Ie(c) 


'(£) Ri'^\e+u) 








;+2u) 


Ru(£+") 




Sf'{E+2u) 
■ to- 


(c) , 




R„(e+u) 




(e+2u) 



Fig. 1.5. (a) Effective circuit representing current conversion at the interfaces of 
the central superconducting island Sc- Resistors, i?p and -Rq stand for an Andreev- 
and a normal processes respectively. The role of voltages at the nodes is played by 
the electron an d hol e distr ibution functions, (b) An illustration of the boundary 
conditions Eqs.(1.34)-(1.36) in terms of a pyramid-circuit is given in this figure. 
Electron and hole currents entering the left side of the pyramid flow in one normal 
layer, the right currents flow in the other normal layer. The effective resistance _Rd 
describes the "direct" quasiparticle transmission from one normal layer to the other 
through the superconductor and the resistance _Rb describes Andreev processes, c) 
Equivalent 3D-sketch of the circuit (b). 



1.7 Spectral current flow through the superconducting 
grains 



The circuit shown in Fig. 1.5 a is the graphic representation of the boundary 
conditions to Eq.(1.2) at the edges of the superconducting island. It is con- 
structed from the circuit units shown in FigT.2a. We consider the case where 
the size of the superconducting island is less than the charge imbalance length, 
and therefore the coordinate dependence of the quasiparticle distribution func- 
tions at the island can be neglected. Solving the Kirchhoff equations for the 



circuit shown in Fig. 1.5 a we exclude the quasiparticle distribution functions 



corresponding to the superconducting island and express the spectral currents 
through the quasiparticle distribution functions in the normal layers: 



nf\e)~n'i\e) n^^\e ^ 2u) ~ n-^\e + 2u) 



.(2) 



Me) 



Ro{e + u) 



i?B(e + u) 



Roie - u) 



R-s{e - u) 



± 



(1.34) 
(1.35) 
(1.36) 



1 Synchronized Andreev Transmission in Chains of SNS Junctions 



13 



The effective resistance describes the "direct" quasiparticle transmission 
from one normal layer to the other through the superconductor and the re- 
sistance -Rb describes the Andreev processes, see Fig. 1.5b. Note that the di- 
rect and indirect transmissions here are different from the so-called "elastic 
co-tunneling" and "crossed Andreev tunneling" |2^, Q processes where Bo- 
goliubov quasiparticles tunnel below the gap through a thin (with the width 
of the order of the Cooper pair size) superconducting layer. The probability of 
these tunneling processes decreases exponentially if the width of the supercon- 
ducting layer exceeds the Cooper pair size, and they occur without generating 
supercurent across a superconductor (the supercurrent flows "virtually"). The 
size of superconucting islands of the SNS arrays that we consider here exceed 
well the Cooper pair size, and the current of the quasiparticles with the ener- 
gies below the hap converts at the NS interface into the supercurrent across 
the S-islands and then transforms again into the quasiparticle current at the 
opposite SN-interface. 



1.8 Recurrent relations 

We have demonstrated that there is a direct correspondence between the effec- 
tive electric circuit and the solution of the Usadel equations with the appropri- 
ately chosen boundary conditions. The effective circuit describing transport in 



SNSNS-array is shown in Fig. 1.6. We choose the direction of the current flow 
in such a way that the electron, /o, and the hole, Ih, currents go in opposite 
directions. The expression for the total current then assumes the form: 

I{V) = -^J de{h + h). (1.37) 

The spectral currents I^. and /h satisfy in general the relation: /o(e) = 
—I\i{e:)\v^-v- Similarly, nc(e) = n\-i{e)\v^-v , ensuring the identity I(—F) — 

Writing down the Kirchofs equations for potential distribution at the cir- 
cuit in Fig. 1.6 we arrive at the recurrent relations, see Appendix |a[ 



7^(£, -u, -y)4(e)™p(°) [e - n)I,{e - 2u)- 

p^^\e)Ue)- p^^\e-V)h{e-V)^n,{e)-n,{e-V), 
7^(e, u, V)h{e) - p(°)(e + u)h{e + 2u)- 

p(^)(£)4(e) - p(^)(e + V)Ue + 2V) = n,{e + V)- n,{e). 

Here the effective resistance TZ = Ri+ R2 + p^^°^^ , where 
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«F(e) 



ni\-2u)i 



/?P™(e-V)/^ 



\£-2V)^ 



-Rj^{£-u)i 



\Rjiie-u) 

/?b(£+m)I 

I n?\e-2u) 

n'^\z+2u)j 



|/?b(£+z<) 7?o(8+m) 



-7?b(£+m\ 2 



4(£-2m) /(^(8+2m) 



7 p(2RV, 



'(e) 



-(2R)/ X ' 



'/h(e+2V) 



«f(S-V) 



«p(£+V) 



Fig. 1.6. MAR in a SNSNS array. The graph shows the efTective circuit for quasi- 
particle currents 1^ and /h in the energy space. The role of voltages here play quasi- 
particle distribution functions. Boxes, triangles and ovals play the role of effective 
resistances that come from Usadel equations and their boundary conditions. 



.,([>o<) 



p = (1/2) E + + CVu + Rt%}. 

a=± 



- R 



il/2){R\ 
p(^) = (l/2){4^"^-i?L'"^}, 
p(^) = (l/2){i?^'^^ 



(1.40) 

(1.41) 
(1.42) 
(1.43) 



In the normal state of the array (or if |e| 3> /\) 7?. reduces to a normal resis- 
tance of the array whereas p^^-* and p^^"^ vanish. Then we find from Eqs.(1.38)- 
( PTSSI ) that 4(e]_= [n^ie) - n^(e -V)]/n, and /c(e) = K(e + F) - nF(e)]/7^ 
that with Eg. (1.37) reproduces the Ohm's law, X = V/TZ. 



It is easy to find the island potential in the case of symmetrical array when 
the transmitivities of the island-normal metal interfaces are equal as well as 
the transmitivities of the lead- normal metal interfaces and Ri = i?2. Then 



the resistances R 



(il) _ p(lC) _ p(2C) 



R'r' = R\ 



and for the symmetry reasons. 



V/2. At the same time the recurrent relations Eqs.(1.38)-(1.39) become 
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invariant under the substitution /e(e — V)— Ihis) and reduce to the relation: 

n{e, V)Ue) - p^^\e)I,{e -V)- p^^^e + V)h{e + V)^ n^e + V)- n^{e), 

(1.44) 

where 

7^(£, V) = 7^(£, y/2, V) - p(°)(e + V/2) = 



R^e) + (1/2) ^ + J, (1.45) 



where R^{e) = i?i + i?2 + R^-':J^+ 



The recurrent relation, Eq.(1.44), is similar to that of a (symmetric) SNS 
junction, see Ref. |lj, ^ and Appendix p|, b ut in our case the normal re- 
sistance Rn{£) becomes energy dependent |22(|. In other words, a symmetric 
SNSNS array has the same transport properties as a single (symmetric!) SNS 
junction, but with the energy dependent resistance of the normal layer. The 
imbalance resistance Ryf^ has singularities at the energy corresponding to 
the gap edges of the superconducting island in the center of our SNSNS ar- 
ray. This is the origin of the subharmonic singularities in the current- voltage 
characteristics at voltages 2A/V = n/2, n = 1, 2, . . ., contrasting the "conven- 
tional values" in an S NS junction determined by the relations 2A/V — n. It 



follows from Eq.(1.45) the unusual subharmonic singularities should disappear 



if the resistance of the normal layer greatly exceeds the resistance of the SN 
interfaces. Then Ri R_ + ^--,e+v/2 ^.nd the central superconducting 
island of the SNSNS array effectively "disappears" and the array completely 
transforms into a SNS junction [p2|. 



1.9 Results and Discussion 

Calculation of the current-voltage c harac te ristic s 1{V) requires numerical 
solving of the recurrent relations, Eqs. ( 1.38|) -( 1.39 ). To accomplish the numer- 



ical task, we have developed a computational scheme allowing to bypass insta- 
bilities caused by the non-analytic behavior of the spectral currents /o(h)(£), 
which poses the major computational challenge. The procedure is as follows: 
first, we fix certain chosen energy e and identify the set of energies connected 
through the equations in the given energy interval, solving afterwards the 
resulting subsystem of equations. We then repeat the procedure, until the 
required energy resolution of 5e = 10~^Z\ is achieved. Typically, up to 10^ 
linear equations had to be solved for every given voltage, but the complexity 
of the coupled subsystem depends on the commensurability of u and V . 



Figure 1.7 shows the comparative results for the SNSNS junction and 
two SNS junctions in series. The latter corresponds to the case where the 
size of the central island well exceeds the energy relaxation length, Lc > ^e- 
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I I I I I I I I I I I I I I I 1 1 I I I I I I I I I I I I I I I I I 




I I I I I I I I I I I I I I I I ' I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — H 

1.5 2.0 2.5 3.0 1.5 2.0 2.5 3.0 
eV /A 

tot 

Fig. 1.7. Left panel: DifTerential resistances as functions of tiie applied voltage Vtot 
(around n = 1 in Eq. (1)) for the SN1SN2S junction. The fractions 3/4, and 4/5, and 
9/10 represent the ratios of resistances of the normal regions, R1/R2. The differen- 
tial resistance dV/dl of the SN1SN2S junction demonstrates the pronounced SAT 
spike at Kot ~ 2A/e, irrespectively to the partial voltage drops. The SAT spike is 
sandwiched between the two additional spikes corresponding to individual MAR pro- 
cesses occurring at junctions SNiS and SN2S for m\,m2 — 2. The voltage positions 
of these features depend on R1/R2. Right panel: The corresponding dV/dI{Vi + V2) 
for the two SNiS and SN2S junctions in series as they would have appeared in the 
absence of the synchronization process, i.e. in the case where Lc > £e- These dV/dl 
dependencies were calculated following (with transmissivity W=l). 



We display the differential resistances as functions of the applied voltage, 
which demonstrate the singularities in Andreev transmission more profoundly 
than the I-V curves. There is a pronounced SAT spike in the dV/dl for 
an SNSNS junction at Vtot = 2A/e. The spike appears irrespectively to the 
partial voltage drops in the normal regions and is absent in the corresponding 
curves representing two individual MAR processes at the junctions SNiS and 
SN2S. 

The resonant voltages of the SAT singularities can be found from the 
consideration of the quasiparticle trajectories in the space-energy diagrams. 
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Such a diagram for the first subharmonic, n — 1 and ratio R1/R2 = 3/4 is 



given in Fig. 1.1. A quasiparticle starts from the left superconducting electrode 
with the energy e = — Z\ to traverse Ni, and the quasiparticle that starts from 
the central island Sc with the same energy as the incident one to take up upon 
the current across the island N2, and hit Sr with the energy e — A (the ABCD 
path, the corresponding path for the hole is D'C'B'A'). In general, relevant 
trajectories yielding resonant voltages of Eq. (1.1) have the following structure: 
they start and end at the BCS quasiparticle density of states singular points 
(e = ±Z\), contain the closed polygonal path, which include MAR staircases 
in the normal parts and over-the-gap transmissions and Andreev reflections, 
and pass the density of states singular points at the central island. Apart 
from the main singularities [Eq. (1.1)], additional SAT satellite spikes appear 
at V ~ {2A/e){p + q)/n, where p/q is the irreducible rational approximation 
of the real number r = R1/R2, (we take i?i < R2), and n ^ {p + q). 

The achieved qualitative understanding enables us to observe that the 
manifestations of the SAT mechanism in an experimental situation becomes 
even more pronounced with the growth of the number of SNS junctions in the 
system. To see this, let us assume that the resistances of the normal islands in 
a chain of SNS junctions are randomly scattered around their average value 
i?o and follow Gaussian statistics with the standard deviation ~ crRo, 
where a is dimensionless. Accordingly, the dispersion of the distribution of 
the MAR resonant voltages is characterized by the same a, and the MAR 
features get smeared. Let us distribute the voltage drop 2A/e among the 
n successive islands. Then the quasiparticle SAT path starts at the lower 
edge of the superconducting gap at island j, traverses n — 1 intermediate 
superconducting islands and hits the edge of the gap at the j + n-th island 
in the chain. The standard deviation of the voltage drop on the n islands 
grows as ^/n resulting in a voltage deviation per one island oc l/\/n, i.e. the 
dispersion of the distribution of Vn drops with increasing n: (Tsat — <y/\/n. In 
contrast to the MAR-induced features, with an increase of n, the subharmonic 
spikes at voltages Vn per junction due to SAT processes become more sharp 
and pronounced. 



1.10 Conclusions 

In conclusion, we have developed a nonequilibrium theory of the charge trans- 
fer across a central superconducting island in an SNSNS array and found that 
this island acts as Andreev retransmitter. We have shown that the nonequi- 
librium transport through an SNSNS array is governed by the synchronized 
Andreev transmission with the correlated conversion processes at the opposite 
NS interfaces of the central island. The constructed theory is a fundamental 
building unit for a general quantitative description of a large array consisting 
of many SNS junctions. 
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A Recurrent relations for the quasiparticle currents and 
distribution functions 



The Kirchhoff's laws applied for the circuit in Fig. 1.6 generate the following 
linear system of equations: 

-(2R)_-(2R) _ T/^ _ 

h{e-2V)^ -t- , (1.46) 

^5(2R) _(4) / 

^h(£) - — ^iW) -I ^iW) ' y^-^') 

-(2) _ (2R) 

hie) = ^ „ " , (1.48) 

_(1) _ _(2) _(1) _ _(2) 

4(e) = + ^ , (1.49) 

_(2) _ _(1) -(2) _ _(1) 

Ic{e-2u)^— — g '■ + ^ , (1.50) 

-(IL) _ -(IC) 

h{e) = "h^^h^ (1.51) 
Hi 



n 



and 
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^(2H) 



5(2R) ^(2R) 



h{e + 2V) 
/h(e + 2u) 



'■h,E+2V 



n^{e + V) 



p(2H) 



p(2R) 



= (2n) 



''h,e+2V "-ce 
p(2R) 

7(2R) _ f)(2C) 

''e,e '^e,e 



p(2R) 



i;(2C)_^(lC) ^(2C) _ f)(lC) 



n(12) 

-(IC) -(2C) 



R 



(12) 



"'li,e+2« "'h,e+2u '^c,s "'c,s 



R 



(12) 



i? 



(12) 



i?l 



5(1L) 



nj,(e) 



p(iL) ' n(iL) 



(1.53) 
(1.54) 
(1.55) 
(1.56) 

(1.57) 

(1.58) 
(1.59) 



Eqs.( 1.4^ )-( |l.59D are the recurrent relations (i.e. the relations coupling the 
functions at energy e with the functions at e ± y) for the currents and the 
distribution functions. 

It follows from Eqs.(|l.48|),(P^) that 



i^e [R + R] = - fir + - 



(1.60) 



The distributions functions ente ri ng Eg . ( 1.60 ) we can express below through 
the currents. Combining Eqs.( 1.49 )-( 1.50| ) we get, 



-(ic) -(2c) \_^h,e-^i^€—u -^c.e — 2u-^D,e~u\ -^d,£— u -^B.e— u 

At the same time from Eqs.(|l.52|),(|T59|) follows that 



(1.61) 



.(IL) 



np(£) + R 



(IL) - 



n(iL) I n(iL) 



(1.62) 



and finally from Eqs.( |l.4^ )-( 1.47 ) we get 



= (2R) _ „ /'e- T/A I n(2r 



(2R) 

e-V 



-io,s-2v Ri'^u + {R+:-i-v + ^: 



•■p.E-y 



r)R(2R) I n(2R) 



(1.63) 



Combining Eq.( |L60|) a nd Eqs.(p^-(|r6^) we find the recurrent relation 
for the currents, Eq.( 1.38 ). Similar procedure helps to derive Eq.( 1.39 ). 
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B Charge transport in SNS junctions 



We discuss below the transport properties of SNS and SNN' junctions to make 
a mapping between our technique and the well-known results obtained before 
us. 



The recurrent relations, Eqs.( 1.3q )-(1.39), solve the transport problem in 
a SNS junction in the incoherent regime. Then there is no island, so p^°'' — 
and we should remove the island resistances with the indices (IC) and (2C) 
from the coefficient functions of the recurrent relations. So, 



(1.64) 



7^(e, V)I,ie) - p^^\e)h{e) - p^^>ie + V)h{e + 2V) ^n^e + V)- n^{e). 

(1.65) 

where, for example, 

TZie, F) = i?i + i?2 + (1/2) J2 {RtJ + ^Siv}- 



(1.66) 



Eqs.(1.64)-( ]l.65 ) are invariant under the following transformation, /c(e - 



V) — >■ /h(e), if at the same time we exchange the resistances, R±^^ O i?''^'^^ 



Thus the relation, /g {e — V) = Ih (e) and the reduction of the recurrent rela- 



tions to the one equation for I,, or for Ih as it was done in Ref. |15 : 

7^(e, V)I,is) ~ p^^\e)h{e - F) - ^(^'(e + V)Ue + F) = npie + F) - ^^(e) , 

(1.67) 



holds only for a symmetric SNS junction with K 



(IL) 



R 



(IR) 

±.e ■ 



To sum mari ze here our consideration [summarized by the recurrent rela- 
tions Eqs.( 1.64 )-( 1.65 )], reduces to that presented in Ref. ||l5[ only in the 
case where the contacts are symmetric and the assumption /o(e — V)— -^h(e) 
holds. 
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